Effect of the s,p-d exchange interaction on the excitons in (Zn,Co)0 epilayers 
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We present a spectroscopic study of Zni-xCoajO layers grown by molecular beam epitaxy on 
sapphire substrates. Zni-iCo^O is commonly considered as a promising candidate for being a 
Diluted Magnetic Semiconductor ferromagnetic at room temperature. We performed magneto- 
optical spectroscopy in the Faraday configuration, by applying a magnetic field up to 11 T, at 
temperatures down to 1.5 K. For very dilute samples (a; < 0.5%), the giant Zeeman splitting of the 
A and B excitons is observed at low temperature. It is proportional to the magnetization of isolated 
Co ions, as calculated using the anisotropy and g-factor deduced from the spectroscopy of the d-d 
transitions. This demonstrates the existence of spin-carrier coupling. Electron-hole exchange within 
the exciton has a strong effect on the giant Zeeman splitting observed on the excitons. From the 
effective spin-exciton coupling, {No{a — P)}x ~ 0.4 eV, we estimate the difference of the exchange 
integrals for free carriers, TVoja — /3j ~ 0.8 eV. The magnetic circular dichroism observed near the 
energy gap was found to be proportional to the paramagnetic magnetization of anisotropic Co ions 
even for higher Co contents. 

PACS numbers: 75.50.Pp, 75.30.Hx, 78.20.Ls, 71.35.Ji 



I. INTRODUCTION 

Theoretical predictions, ^•^•''•^ and first reports of room 
temperature ferromagnetism in Zni_2;Co2;0)SiSiLS424i2iii 
have strongly stimulated researches on this wide bandgap 
Diluted Magnetic Semiconductor (DMS). However, the 
origin of the observed ferromagnetic behavior remains 
controversial, and the existence of a coupling between the 
magnetic and electronic properties has still to be demon- 
strated, and its magnitude evaluated.^ 

Most of reported studies rely on magnetometry. Direct 
magnetization measurements on thin DMS layers are par- 
ticularly difficult, because the magnetic signal from the 
thin layer is usually small compared to the diamagnetic 
contribution from the much thicker substrate. Magneto- 
optical spectroscopy on the excitonic transition and Mag- 
netic Circular Dichroism (MCD) in the vicinity of the 
bandgap can be used very efficiently to measure the tem- 
perature and magnetic field dependence of the magne- 
tization, using the giant Zeeman effect of substitutional 
magnetic ions coupled to the band electrons, and to dis- 
criminate the contribution from magnetic ions incorpo- 
rated in spurious phases. In addition, the observation 
of the intra-ionic d-d transitions gives us an information 
on the electronic structure of the incorporated magnetic 
impurity, and allows us to measure accurately the param- 
eters describing its ground level and governing the mag- 
netization of the isolated impurities. Finally, combining 
the two sets of data, on can deduce the magnitude of the 
s,p-d exchange interactions between localized spins and 
free carriers. 

This paper describes the results of such a study con- 



ducted on a set of Zni-ajCo^^O layers with various Co 
concentrations, including very dilute samples which fea- 
ture well resolved excitonic lines, so that we can measure 
the exchange Zeeman splitting. It fills the gap existing 
between recent magneto-optical studies of Zni_j;Co2;0 
layers with a high Co content, where large linewidths 
prevent any direct observation of the exchange Zeeman 
splittingjiiiiSiiiiiiiSiiSiii and early studies of very dilute 
Zni_a;Co3;0 bulk samples without applied field liSiiSiSfliSi 



II. SAMPLES AND EXPERIMENT 

About 1/im-thick layers were grown on sapphire sub- 
strates by plasma-assisted molecular beam epitaxy. Two 
dimensional growth was achieved for a growth tempera- 
ture of 560°C (i.e., 50°C higher than the optimal growth 
temperature used for ZnO), resulting in streaky RHEED 
patterns. The Co content of several layers was measured 
by energy dispersive x-ray analysis (EDX). For low Co 
contents {x < 1%), the full widths at half maximum of 
the x-ray rocking curves are in the range of uj ~ 0.15° 
along (002), (105), and (105). The low, identical values 
of ui measured both for (105) and (105), indicate a large 
column diameter, close to 1/im (up to a; = 15%). While 
the column diameters remain large, uj values are found 
to increase slightly and gradually with the Co concentra- 
tion. For X = 15%, measured uj values are 0.32°, 0.22°, 
and 0.66° along (002), (T05), and (101) respectively. For 
all compositions, the c-axis of the wurtzite structure is 
perpendicular to the film plane. No other orientation 
or column rotation were detected. The conductivity of 
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the films is n-type, with residual carrier concentrations 
He < I.IO^^ cm~^, below the Mott transition. The thick- 
ness of the layer was checked on the electron microscope 
image of a cleavage plane. 

Reflectivity, transmission and photoluminescence (PL) 
measurements were performed in the Faraday configura- 
tion, with the optical axis and a magnetic field (up to 
11 T) both parallel to the c-axis, at temperatures down 
to 1.5 K. A high-pressure Xe lamp was used for trans- 
mission and reflectivity; PL was excited with a He-Cd 
laser. 



III. SPECTROSCOPY OF COBALT IONS 

We measured the characteristic absorption lines and 
bands in every studied Zni_j;Coa;0 sample from the most 
diluted one (x — 0.1%) up to the most concentrated one 
{x = 35%). Examples are given in Fig. n(a) for the 
spectral range of interest in the following. These lines 
and bands have been already identified as intra-ionic d-d 
transitions of cobalt in bulk samples by P. Koidl»ii For 
substitutional Co ions, the different eigenstates associ- 
ated to the 3d^ configuration are labeled using the nota- 
tion of McfarlaneiSS 

We will particularly study the transitions between the 
ground state and the "^E excited state arising from 
the cubic part of the crystal field. Both quadruplets 
are further split into two components (two Kramers dou- 
blets) , 2 A and E, by the spin-orbit coupling and the trig- 
onal component of the crystal field. For the ground 
states, we label |±i) the two spin sublevels of E and 
|±|) the two spin sublevels of 2A. The 'E - 2A sphtting 
of the ground states is to small to be resolved in the 
spectra of Fig. ^a), while that of the ^E excited state is 
clearly visible. 

Other lines which can be seen at higher energy in 
Fig. n^a) have been attributed by Koidl^® to transitions 
to other excited levels. They are broader, and hence more 
difficult to study quantitatively in thin layers since trans- 
mission spectra are plagued by interferences between the 
surface of the layer and the interface with the substrate. 

In this section we use the '^A2 ^ "^E transition to (A) 
estimate the Co content in samples which have not been 
characterized by EDX (particularly at low Co content), 
(B) confirm the parameters describing the ground state 
and the excited state (anisotropy and g- factors), and the 
oscillator strengths of transitions and (C) from a com- 
plete description of the intensities of the absorption lines, 
confirm the population distribution in the ground level 
and deduce the magnetization of the system of localized 
spins. 

A. Calibration of Co concentration 

In a set of samples, where the cobalt concentration 
was determined by EDX, the absorption coefficient of 
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FIG. 1: (a) Absorption spectra of Co^^ in Zni_a;Coa;0 sam- 
ples with various cobalt concentration x, at T = 1.6 K. Two 
zero-phonon intra-ionic transitions are identified: *^A2 — > 
'^EE (at 1875 meV) and '^A2 '^E2A 1880 meV). 
The absorption coefficient of the — > ^E2A transition is 
shown by a dotted horizontal line pointing to the right figure, 
(b) The absorption coefficient at 1880 meV, as a function of 
the Co concentration determined by EDX (full circles). The 
observed linear dependence allows us to determine also the 
Co concentration in very diluted samples (open squares and 
vertical arrows), where the EDX technique is not sensitive 
enough. 



the ■^A2 '^E2A fine [at 1880meV, Fig. [Ha)] increases 
linearly with the Co concentration up to 6%. [full circles 
in Fig. n^b)] . The calibration curve obtained with these 
samples was used to determine the cobalt concentration 
in other samples [open squares in Fig. ^b)] . 

The linear behavior observed up to a Co concentration 
of about 6% indicates that the absorption coefficient is 
sensitive to the total concentration of Co^+ ions substi- 
tuting Zn in the wurtzite lattice. In particular, it does 
not discriminate the different nearest neighbor configu- 
rations: isolated ions, antiferromagnetic or ferromagnetic 
cobalt pairs, triplets, or other complexes. In the cation 
sublattice of the wurtzite structure, every cation has 12 
nearest neighbors: 6 neighbors in the same layer, and 3 
neighbors in the next layer on each side, at almost the 
same distance. These 12 neighbors are usually considered 
as equivalent (see, e.g., the study of Zni_a;Mna;0^). As- 
suming a random distribution of Co ions, the probability 
to find an isolated cobalt ion (without any cobalt ion in 
a nearest neighbor position) is given by (1 — x)^^. For a 
total concentration x — 5.6%, half of the Co ions have at 
least one Co ion as a nearest neighbor. That means that 
the concentration of isolated Co ions, Xs — x{l — x)^^, 
significantly deviates from the straight line drawn in 

Fig.nb)]. 

The magnetic properties of a set of isolated Co ions and 
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(a) Photon Energy [meV] (b) Photon Energy [meV] (d) Magnetic Field [T] 

FIG. 2: Magneto-optical spectroscopy of tlie *yl2 ^ ^EE transition in a Zno.98Coo.02O at T = 7 K, in a'^ (top) and a~ 
(bottom) circular polarizations. The magnetic field up to 11 T is parallel to the c-axis and to the propagation of light (Faraday 
configuration), (a) PL spectra, (b) Absorption spectra, (c) Position in energy of the observed PL (full circles) and absorption 
(triangles) lines, and (solid lines) values calculated using the parameters of Table ITTl (d) Energy level diagram, as a function of 
the applied field, and transitions. The spin levels of the ground level quadruplet are marked. 



of nearest neighbors complexes are expected to be very 
different at low temperature. Antiferromagnetic nearest- 
neighbor pairs are blocked antiparallcl; in many DMSs, 
this gives rise to a series of steps and plateaux in the 
high field part of the magnetization cycles. At smaller 
fields, the magnetization is well described if one replaces 
the total number of spins, x, by an effective number of 
spins, Xeff- The value of x^ff has been determined in 
a large number of DMSs with the zinc-blende structure. 
It does not depend on the chemical nature of the DMS, 
which suggests that it corresponds to a statistical distri- 
bution in an ensemble of randomly occupied sites. For 
small values of x, the above value of Xs coincides with 
the experimental value24 of Xeff in Cdi_2:Mna;Te. We 
shall use this approximation in the following. 

B. Magnetospectroscopy 

We combine PL and absorption data involving the 
<-> "^EE intra-ionic transitions in order to check the 
parameters governing the evolution under magnetic field 
of the '^A2 ground state. Fig. Ol^a) shows the spectra 
around the zero phonon PL lines near 1875 meV. Two 
PL lines with equal intensities are observed in zero mag- 
netic field. They have been identifiedi^ as transitions 
from the excited state "^EE to the two components of 
the ground state (|±^) and |±|)). In the presence of a 
magnetic field parallel to the c-axis, we observe six PL 
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TABLE L Optical selection rules for dipole transitions in trig- 
onal symmetry (see Macfarlane^S). These selection rules di- 
rectly reflect the conservation of total momentum, modulo 
3. The TT polarization corresponds to the active field of the 
light parallel to the c-axis of the wurtzite crystal and can- 
not be observed when light propagates along the c-axis; the 
(7 polarization corresponds to the active field perpendicular 
to the c-a^xis, and is observed with both helicities when light 
propagates along c. 

lines [see Fig. El^a)], which were identified as transitions 
between the different spin sublevels of the ground and ex- 
cited states split by the Zeeman effect. In agreement with 
the selection rules for dipolar transitions (see Table P))^ 
for light propagating along the c-axis, two transitions are 
forbidden, three lines are observed in polarization and 
three in <t~ polarization. 

The transitions observed in PL are also visible in ab- 
sorption [Fig. I3b)]. Fig. |3c) shows the position of the 
PL and absorption lines as a function of the intensity 
of the magnetic field up to IIT. The solid lines are cal- 
culated using 4 fitting parameters: values of the Lande 
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Parameter 


This work 


References 


'^EE 


9\\ 


-3.52 


-3.358'" 




9\\ 


2.28 


2.2384* 


l±|)-l±|) 


2D 


0.69 meV 


0.68 meV* 


|±|) ^ ^EE 


energy 


1875 meV 


1877 meV^ 



"Theory, calculated by R. M. Macfarlane^^ 

'EPR experiment done by N. Jedrecy et aliSI., diluted sample 

■^Optical absorption measured by P. Koid^£, diluted sample 



TABLE II: Parameters describing the spin splitting of the 
ground state quadruplet and the '^A2 <-> '^EE transition. The 
values were derived from PL and absorption measurements of 
Zni_:,Co:,0 with 2% Co. 

factors 5|| — 2.28 for the ground state and g'j — —3.52 
for the excited state, a zero-field splitting of the ground 
state 2D = 0.69 meV, and an energy of the zero field 
transition between |±|) and '^EE equal to 1875 meV. 
These parameters, obtained for a (strained) layer with 
2% Co, are in good agreement with the values obtained 
theoretically or from other experimental techniques (see 
Table InJ . The corresponding energy diagram and the 
scheme of transitions are shown in Fig. |2Id). 

Since we shall be interested in the magnetic properties 
of the Co system, it is important to note that in zero 
magnetic field, the |±i) states are at lower energy than 
the |±|) states. This ordering induces a spin anisotropy 
with an easy axis perpendicular to the c-axis. Due to the 
three times larger Zeeman splitting^^ of the |±|) states, 
if the field is along the c-axis, we observe a crossing of 
|-|) with at 5.2 T [See Fig.|2i;d)]. 
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(a) Energy [meV] (c) Magnetic Field [T] 

FIG. 3: (a) Absorption spectra of Zno.955Coo.045O measured 
in the Faraday configuration, at 1.7 K, in a magnetic field up 
to 11 T, in and a~ circular polarizations. The four tran- 
sitions observed are labeled according to Fig. I^d). (b) Ab- 
sorption intensity measured at T = 1.7 K (symbols) and cal- 
culated using Maxwell-Boltzmann statistics (lines), (c) Mean 
spin of Co calculated from the absorption intensities (sym- 
bols) or from Maxwell-Boltzmann statistics and parameters 
of Table ini (lines). 



C. Populations of the spin sublevels and 
magnetization 

The magnetization of the Co system is determined 
by the population of the different spin sublevels of the 
ground state. We now check these populations from the 
absorption intensities and calculate the expected magne- 
tization. 

For a better accuracy, we used magneto-absorption 
spectra from a sample with a higher Co content, 4.5% 
[Fig- EI 9')] • As each absorption intensity is proportional 
to the population of the initial state, we observe only four 
— > ^EE transitions, which are labeled as in Fig. |21 

Using the parameters determined from the previous PL 
and absorption data (Table ^ , and assuming gaussian 
line shapes with a linewidth equal to 1.44 meV, we were 
able to fit the whole set of spectra with the intensity of 
each line as adjustable parameters. The resulting inten- 
sities are shown by symbols in Fig. E^b). At low temper- 
ature (T = 1.7 K), two transitions (L2+ and L3— ) keep 
zero intensity for every magnetic field. 

The integrated intensity of each absorption line is pro- 
portional to the occupancy of the initial state and to 



the oscillator strength of the transition. Relative oscil- 
lator strengths are easily determined from extreme cases 
where one level is populated. At T = 1.7 K and i? = 11 T 
only transitions L2— and L3-I- are visible [Fig. |31 (a) and 
(b)] because only the |-|) state is occupied. Similarly, 
in zero magnetic field only the and states are 
occupied so that only LI -I- and LI— are observed. Mea- 
suring the relative intensities in zero magnetic field and 
with S = 11 T leads to the conclusion that the relative 
oscillator strengths of L1-I-, L2— , L3-I- are 1 : 0.7 : 0.35 
[see Fig. Elb)]. We expect that the pairs of lines (Ll-I- 
and LI—), (L2+ and L2— ), and (L3+ and L3— ), have 
equal oscillator strengths because of time reversal sym- 
metry in zero field. Independently, the relative oscillator 
strength of L1+, L2— and L3— can be estimated from 
the intensity of PL lines in i? = 11 T, when only the 
lower component {E+) of "^EE excited state is occupied 
during the lifetime [Fig. [^a)]. The oscillator strengths 
are summarized in Table ITlll 

The occupancy of the different spin sublevels can be 
calculated using the parameters of Table ^ and assum- 
ing a Maxwell-Boltzmann distribution. The solid lines 
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Transition 
Excited state 
Ground state 



L1+ LI- L2+ L2- L3+ L3- 

E— E— E— E-\- 

i-i) i+i) Ki) i-i) H) i-i) 



Relative OS (from abs.) 1 
Relative OS (from PL) 1 



1 0.7" 0.7 O.SS*" 

0.8 O.IS*" 



"from time reversal symmetry 
''not very accurate, see Fig. I^b*) 



TABLE IIL Ratio of the oscillator strengths (OS) of the spin 



split components of the A2 



E2A transition in a polar- 



ization. OS equal to 1 corresponds to a peak intensity of 
0.14//^m with a linewidth of 1.44 meV measured on a sample 
with 4.5% Co. 



in Fig. |3Ib) show the corresponding intensity calculated 
for each line using the relative oscillator strengths of Ta- 
ble Experimental data for L1-I-, LI—, and L2— are 
in good agreement with the calculation. A significant 
discrepancy appears for L3-I-. For this line, the determi- 
nation of the integrated intensity is made difficult by the 
overlap with L1-I-. Moreover a magnetic field dependence 
of the oscillator strength cannot be ruled out experimen- 
tally 

The experimental absorption intensities divided by the 
relative oscillator strengths give us a direct information 
about the expected cobalt mean spin, which can be cal- 
culated using 



z) spectro 



\{l2--h+)/r2,i 



1+ 



+ h+)lr2, 



(1) 



where /jv denotes the experimental intensity of absorp- 
tion line Ltv, and r2_i denotes the oscillator strength ratio 
of lines L2± and Ll± (it is equal to 0.7, see Table HiH . 
We do not use L3± since its experimental intensity is 
doubtful. The expected value is shown by symbols in 
Fig. OIc), as a function of the magnetic field along the 
c-axis, for two temperatures: 1.7 K and 6 K. It clearly 
deviates from an isotropic Brillouin function. 

Here again, assuming a Maxwell-Boltzmann distribu- 
tion and using the parameters of Table ITll we can calcu- 
late the mean spin of an isolated cobalt, as a function of 
temperature T and magnetic field along the c-axis. 



{-s^ = 



\ sinh(i(5) -f |e sinh(|(5) 
cosh(i(5) + cosh(|(5) 



knT 



(2a) 



(2b) 



where ks denotes the Boltzmann constant and /is the 
Bohr magneton. Parameters g\\ and 2D are given in Ta- 
ble |n] The result is shown as solid lines in Fig. ^c) . 

Particularly for very low temperatures, the field depen- 
dence of the magnetic moment deviates from the Bril- 
louin function B3/2 characteristic for an isotropic spin. 



A step is visible at 5 T, as an effect of the crossing be- 
tween |-|) and [Fig. Old)]. Also, if we plot the 
inverse of the low- field susceptibility as a function of tem- 
perature, we observe a clear deviation from a Curie law. — 
These two effects mimic the magnetization steps and 
Curie- Weiss law observed to be due to spin-spin inter- 
actions in the case of DMSs containing the isotropic Mn 
spin; here they are simply due to the single-ion anisotropy 
of the Co spin. 

The magnetization (Mz) of the cobalt system is di- 
rectly related to the previous mean spin value [Eqs. (Q) 
or ©]: 



(Mz) = ~g\\ ^J.B Nq Xeff (Sz) 



(3) 



In section HVl we describe results of magnetooptical mea- 
surements performed at photon energies close to the en- 
ergy gap of Zni_a;Coa;0. In section^ we compare them 
with the magnetization above. 



IV. EXCITONIC SPECTROSCOPY 

In this section, we first recall the main characteristics 
of the spectroscopy of ZnO, a wide gap semiconductor 
with the wurtzite structure and a small spin-orbit cou- 
pling. Then we describe the giant Zeeman effect observed 
on sharp excitonic transitions in very dilute Zni_a;Coa;0 
layers. Finally we turn to our observations of near-gap 
magnetic circular dichroism in layers with a larger cobalt 
content. 



A. Zero field results 

ZnO and Zni_j,Coa;0 naturally crystallize in the 
wurtzite structure. In this structure, the trigonal com- 
ponent of the crystal field and the spin-orbit coupling 
split the valence band in the center of Brillouin zone into 
three valence band edges. These splittings are small and, 
in epitaxial layers, they depend on the residual strain. 
Actually, the position and the symmetry of the valence 
edges are still a mater of controversy even in bulk ZnO. 

Three excitons can be detected in ZnO. They are la- 
beled {A, B, C), in the order of their position in the spec- 
tra, independently of their symmetry. Exciton A, which 
appears at the lowest energy, is associated to the valence 
band edge at the highest energy, and C to the edge at 
the lowest energyi22i2Si22i22i^ Both the A and B excitons 
are observed by reflectivity at normal incidence (fc || c), 
hence in a polarization (electric field of the light perpen- 
dicular to the c-axis of the crystal, i? _L c), while exciton 
C can be observed only in tt polarization (E \\ c). 

These three valence band edges result from the com- 
bined effect of the trigonal crystal field (described by a 
parameter noted Ai) and the anisotropic spin-orbit cou- 
pling (two parameters, A2 = aJ^q/S, A3 = Ag^/S). In 
ZnO, the spin-orbit coupling is much smaller than the 



6 




(a) 3350 
3460 



3400 

Photon Energy [meV] 



3450 




(b) 



1.0 2.0 3.0 4.0 
Cobalt Concentration [%] 



5.0 



FIG. 4: (a) Reflectivity spectra measured at T = 1.6 K with 
incidence angle of 45 degrees. Topmost spectrum is for ZnO, 
other spectra for Zni_a;Coa;0 with increasing Co concentra- 
tion (0.1%, 0.4%, 0.5%, 2%). Labels A, B and C identify the 
three excitons which are visible in ZnO and at low Co content; 
Perot-Fabry oscillations are also observed at low photon en- 
ergy, (b) Spectral positions of the excitonic transition in ZnO 
and TiXi^-xCoxO. Symbols marked as "mixed j4-|-i?" means 
that in the case of samples with a Co concentration higher 
than 2% we cannot resolve excitons A and B, and we plot 
the position of a broad structure. The dashed line shows the 
linear fit for the position of exciton yl as a function of Co 
content. 



trigonal crystal field (more precisely, A3 is much smaller 
than Ai — A2). Hence the trigonal field splits the p- 
likc states which form the top of the valence band into 
a doublet (Fs, in-plane p-states) and a singlet (Fi, out- 
of-plane p-states). As a result of the spin-orbit coupling, 
the orbital doublet is split into two (orbit -|-spin) dou- 
blets (Fyj-s), Fg(5)), and the orbital singlet forms a ^7(1) 
doublet with a very small admixture of the '^7(5) states. 

Exciton C, which is active in tt polarization, is unam- 
biguously ascribed to the F7(i) doublet. The attribution 
of excitons A and B is not straightforward. The splitting 
energies are given byi^ 



Ev 



.-Er 



Ai -H3A2 



'Ai - A2 



7(5) 



'Er, 




2Ai (5) 



2AI (4) 



Lambrecht and co-worker^ have obtained the fol- 
lowing parameters for pure ZnO: Ai = 38 meV, 
IA2I = 4.53 meV and A3 = -3.05 meV. We cite only 
the absolute value of A2, because we do not want 
to suggest any particular energy order for the valence 
states: a positive sign of A2 implies that the valence 
state A has the Fg(5) symmetry(2L2S422i2£ a negative 
sign that A is ^7(5) Mt^^t^ Actually, taking into account 
that the spin-orbit A3 parameter is much smaller than 
(Ai — A2), we obtain the approximate splitting energies: 
^r„„ - i?r,,,, = 2A2, £;r,„) - i?r„i, = Ai + A2, and 

^^7(5) - E^Hl) = Al - A2. 

In thin epilayers with the c-axis perpendicular to the 
plane, reflectivity measurements in tt polarization {E \\ c, 
/c _L c) are particularly difficult. In order to observe exci- 
ton C we performed refiectivity measurements with the 
sample tilted by an angle of 45 degrees with respect to 
the optical axis. Then the selection rule is relaxed, and 
we can observe all three excitons [Fig.^Ja)]. Spectra are 
analyzed using the polariton formalism described in the 
Appendix. 

We observe a systematic increase of the excitonic tran- 
sition energies with the Co concentration in all studied 
samples [Fig.0Jb)]. This indicates an increase of the en- 
ergy gap of T^W-i-xGoxO with x. The energy of exciton 
A is Ea{x) = 3372.6 meV + xx 1690 meV. The energy 
differences between the excitonic lines in zero field are 
Eb — Ea = 7 meV and Ec — Eb = 47 meV, whatever 
the Co content. 

Using the approximate splitting energies above, we ob- 
tain Ai = 50 meV and IA2I = 3.5 meV. The difference in 
Ai between the present layers and values previously re- 
ported for bulk samples'^^ may be due epitaxial strain, 
but also for a good part to the fact that polaritons were 
not taken into account. 

Finally, as we observe excitons, we should take into 
account the electron-hole exchange, which is by no way 
negligible in ZnO. As the value of the parameters de- 
scribing the excitons in ZnO is still a matter of debate, 
we keep for a while this simple description, and we will 
come back to this important point in sectionlVlB. 



B. Zeeman effect 

We performed magneto-optical spectroscopy in the 
Faraday configuration, by applying a magnetic field up 
to 6 T parallel to both the c-axis and light propagation. 
In this configuration, excitons A and B are observed in 
both and a~ circular polarizations. The reflectiv- 
ity spectra in each circular polarization were analyzed 
by taking into account the formation of excitonic polari- 
tons. Polaritons are particulary important in the case 
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FIG. 5: (color online) (a) Reflectivity spectra of ZnO at 
B = T and at i? = 6 T in cr"^ and cr~ circular polariza- 
tions (symbols). The positions of excitons A and B, as de- 
termined from the flt of the _B = T spectrum (solid line), 
are marked by arrows. The Zeeman effect induces a small, 
opposite shift of the two excitons, which changes their over- 
lap and induces a change of the reflected intensity, (b) The 
position of excitons A and B in ZnO, as determined from a 
fits of the reflectivity spectra shown in (a), (c) Reflectivity 
spectra of Zno.9g6Coo.004O at _B = T, and B = 6 T in 
and in a~ polarizations. Solid lines are the fits to the spectra, 
(d) The position of excitons A and B, versus magnetic field. 
At _B = 6 T, the value of the splitting is -1.8 meV for exciton 
A and 1.6 meV for B. 



of ZnO because of the strong coupling between excitons 
and photons. The dielectric function and the reflectivity 
spectra near excitonic resonances have been described in 
detail by Hopfield and ThomasM and by Lagoisi^Si^ The 
application of this theory to our case is described in Ap- 
pendix. 

The transition energies of both A and B and the cor- 
responding values of the Zeeman splitting [Fig. |S] where 
obtained in a two-step procedure. First, the whole set of 
parameters (polarizability aoA,B, position lua,b, width 
Ta,b, the dead layer thickness d, and the non resonant 
absorption contribution ie') were determined by fitting 
the zero field spectra. Then, the spectra under magnetic 
field were fitted by adjusting only five parameters: uja,b, 



Ta,b, and the relative polarizability of excitons A and B 
(keeping constant the sum aoA + aos)- 

We did not observe any excitonic diamagnetic shift, 
as expected due to the small excitonic Bohr radius 
{Rx — ISA) in ZnO^^ (the diamagnetic shift is smaller 
than 0.1 meV at B — 6 T, so it is smaller than 
0.003 meV/T^). 

For ZnO, the Zeeman splitting of excitons A and B 
at 6T remains very small (about 0.1 meV). Such a value 
corresponds to an effective excitonic Lande factor less 
than 0.3. Such a small value of the Lande factors for the 
allowed excitonic transitions in ZnO has been already re- 
ported in the Voigt configuration by Blattner et al.^ 
and more recently by Reynolds et alm^ In this work, the 
experiments are performed in the Faraday configuration 
and they confirm unambiguously the previous observa- 
tions. 

In samples doped with cobalt ions, we observe an en- 
hancement of the excitonic Zeeman splitting. The Zee- 
man splitting of excitons A and B is almost opposite. 
It increases with the Co concentration [Fig. and de- 
creases with temperature. Fig. shows the Zeeman split- 
ting of excitons A and i? as a function of the magnetic 
field up to 11 T, at three temperatures: 1.7 K, 7 K and 
20 K. A saturation appears at the lowest temperature 
and the highest magnetic field. The Zeeman splitting is 
proportional to the mean spin projection calculated for 
isolated Co in section HTll [Fas. 1^1 and (0)]. 

This enhancement of the Zeeman splitting which in- 
creases with the cobalt content, and its proportionality to 
the Co magnetization, are essential for this study. They 
are analyzed in section El as a consequence of s,p-d in- 
teractions. 



C. Magnetic Circular Dichroism 

In samples with Co concentration above 2%, the broad- 
ening of the excitonic lines prevents a direct observation 
of the Zeeman splitting of excitons A and B. In such a 
case, we measured the degree of circular polarization of 
the transmitted fight, — Ia-)/{Ia+ + For the 

sake of comparison with published datajiiiiSiiSii^iiSiiSiii 



we plot the MCD defined as: 



MCD 



180° 
47r 



(6) 



where and are the optical absorption coefficients 
in (T^ and polarizations. Assuming a weak absorp- 
tion and neglecting multiple reflections, the MCD can be 
expressed as 



MCD 



180° 



2tiI /^h 



(7) 



where / is the thickness of the sample, and and /o— 
the intensities of transmitted light in cr~ and a'^ polar- 
izations. The MCD depends on the photon energy and 
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FIG. 6: Zeeman splitting of exciton A, measured at _B = 6 T 
and T — 1.6 K, so that the mean spin of isolated Co is 
(Sz) = 1.2 [Eq. J^J]. The horizontal scale (bottom scale) is 
the absorption coefficient of the — > ^E2A transition of 
Co (see Fig. 0. The top scale (nonlinear) displays the den- 
sity of free spins x^ff expected for a random distribution of 
Co on the Zn sublattice. The solid line shows the giant Zee- 
man splitting {No{a — I3))x x x^jj x (-Sz) calculated with an 
effective (iVo(a - P))x = 0.4 eV. 
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FIG. 8: (a) Magnetic Circular Dichroism (MCD) obtained in 
transmission on four samples with various Co concentrations. 
The MCD signal contains three contributions, one (pointed 
by arrow, most intense in the sample with the lowest Co con- 
centration, 2% Co) close to the bandgap at 3.4 eV, the second 
below the gap,^^ the last one (most intense in the sample with 
the highest Co concentration, 15% Co) in coincidence with the 
internal transitions of Co (1.8 eV - 2.3 eV). Spectra are shifted 
for clarity. 
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FIG. 7: (color online) Dots and left axis: Zeeman splitting 
of excitons A and B measured at three temperatures; lines, 
right axis: mean spin of isolated Co ions as given by Eq. J^J. 
The sample has x^jj — 0.36%. The left and right scales are 
chosen so that {No{a — l3))x ~ 0.40 eV for exciton A and 
0.39 eV for B. 



it is usually very strong near a resonant line or band 
split by Zeeman effect. In simple cases, the Zeeman 
splitting AE can be estimatedi^ from the expression 
MCD = 180747r x dk/dE x AE, where dk/dE is the 
derivative of the absorption coclRcicnt in zero field. How- 
ever, to use properly this equation, the line shape should 
remain constant: this is not the case in Zni_a;Coi:0, 
where the line shape changes significantly, and where the 
A and B excitonic transitions are close together. Addi- 
tionally, they split in opposite directions, as clearly ob- 
served on the reflectivity spectra of Fig.|Sl These changes 
of the lineshape are the reason, why we did not attempt 
to determine the absolute value of the Zeeman splitting 
from these MCD measurements. In each sample, we only 
analyzed the temperature, spectral and magnetic field 
dependence of the MCD signal. 

We observe three main MCD features over a large spec- 
tral range, from 1.7 eV to 4 eV (Fig. 

The first one is related to the energy gap (3.4 eV). In 
samples up to a Co content of 11%, it is seen as a sharp 
negative component, which broadens as the Co content 
increases. It is accompanied by a broader positive com- 
ponent from 2.9 to 3.2 meV, which was attributedii*iSiii 
to a charge transfer transition, or level to band transi- 
tion. It is this peak which was used to study the onset 
of fcrromagnetism upon annealing^Siii 

The third one corresponds to the internal transitions 
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FIG. 9: (color online) Symbols, left axis: Magnetic Circular 
Dichroism in transmission, at a photon energy 3395 meV, 
close to the bandgap; lines, right vertical axis: mean spin of 
isolated Co ions [Eq. @]. The sample is Zni_2,Co2,0 with 
2% Co. 

of cobalt, and it is observed between 1.8 eV and 2.3 eV 
in samples with a Co concentration over 0.4%. Its inten- 
sity systematically increases with the Co concentration. 
This MCD signal is quite complex and mainly results 
from the different intra-ionic absorption lines analyzed 
in section IITTI 

The temperature and magnetic field dependence of the 
MCD signal near the energy gap was found to be propor- 
tional to that of the magnetic moment of isolated Co ions 
[Eq. 10)], even in a sample with a cobalt content of 4.5% 
[Fig. 13 . For higher concentrations, the resonant signal 
becomes too weak to be detected with reasonable accu- 
racy by our experimental setup. 

V. ANALYSIS 

The magneto-optical effects observed near the energy 
gap and described in sections IIVI B and II VI C are now 
interpreted as a consequence of the s,p-d interaction 
between Co ions and band carriers. We first quantify 
the ion-exciton coupling (subsection A); then (subsection 
B) we discuss the effect of the electron-hole exchange in 
the exciton; finally (subsection C), we discuss briefly the 
magnetic properties of the present Zni_2;Coa;0 layers, as 
they appear from magneto-optical spectroscopy. 

A. Excitonic Giant Zeeman splitting 
without electron-hole exchange 

In II- VI DMSs, the presence of the s,p-d exchange in- 
teractions leads to large magneto-optical effects (giant 



Zeeman splitting and giant Faraday rotation). These 
effects have been studied in detail in several wurtzite 
II-VI DMSs like Cdi_^Mn^Se,Si^ Cdi_^Co^Se,iMi 
Cdi_^Fe:,Se,'*2 Zm.^^Mn^Se^^ and Cdi_^Cr^S4i 

The effective Hamiltonian describing the giant Zeeman 
effect is 

Ff/ = -iVoaXe//(S).s, (8) 
for the conduction band, and 

</ = -iVo/?Xe//(S)-S, (9) 

for the valence band. In these expressions, iVo is the num- 
ber of cations per unit volume, Xeff the free magnetic 
ion content, a the exchange constant for the conduction 
band, (3 the exchange constant for the valence band, (S) 
the mean spin of the free magnetic ions, in our case Co 
ions, and § is the spin of the carrier. The mean value of 
the isolated Co spin (Sz) has been already obtained as a 
function of temperature in section III [Eqs. JQ), ©I, for 
a magnetic field applied parallel to the c-axis. 

In order to calculate the effect of these effective Hamil- 
tonians, we must take into account the structure of 
the valence band. We note Ge = l/2NQaXeff{-Sz) and 
Gh = l/2No(3xeff{-Sz)- Neglecting the contribution of 
the weak Zeeman effect of the host semiconductor, the 
energies of the different excitonic transitions under mag- 
netic field are given by^SmiLiSiiMi 

E^f^^^=Eo-Ai-A2TGe±Gh, (10) 



i?r-i, =£^o-^^y^±Ge-i5±, (11a) 



E± = ^ (^^Y^ ± 4- 2Ai, (lib) 

where {Eq — Ai — A2) is the zero field energy of the Fgjsj 
excitonic transition. We do not observe any feature re- 
lated to the C exciton (r7(]^)) in the configuration k \\ c, 
as a consequence of the small spin-orbit coupling in ZnO 
[A3 << (Ai — A2), see section Hvl A]. Hence we can sim- 
plify Eq. dm to E^^ = Eo - Ai + A2±GeT Gh. 

7(5) 

Such a simplification leads to opposite Zeeman splittings 
for the two optically active excitons : E^'^ — E^~^ = 

-{E^+^^ - E^-^J = 2(Ge - Gh) = No{a - (3)x,ff{-Sz). 
This is what we observe experimentally. One conse- 
quence however is that the a and /3 exchange integrals 
cannot be determined from the experiment: only their 
difference can be obtained. Another consequence is that 
the contributions of the two active excitons to the MCD 
tend to cancel each other in samples where the linewidth 
is not small with respect to the A-B splitting. 

There is also a direct relationship between the sign of 
the giant Zeeman effect and the valence band ordering. 



10 



From the sign of the giant Zeeman spHtting we observe 
on A and B, a conventional valence band order (Fgj^s), 
^7(5), r7(i)) implies a positive sign of a — /? (as usually 
observed in II- VI DMSs). This ambiguity would disap- 
pear in more concentrated samples, where the excitonic 
Zeeman splitting will become comparable to the differ- 
ence (Ai — A2), if the lines remained sharp enough. 

In our diluted samples {x < 0.5%), the magnitude of 
the exchange integral difference determined on the op- 
tically active excitons, {NQ{a — l3))x, is about 0.4 eV 
[Fig. El . 

In samples where the two excitons are not resolved 
(such as in our samples with a few % Co), one can at- 
tempt to analyze the MCD phenomenologically, by de- 
ducing a Zeeman splitting from a comparison between 
the MCD spectra and the derivative of the zero-field ab- 
sorption spectrum. We obtain this way (not shown) a 
value of the coupling which is of similar magnitude, or 
slightly smaller. As noticed above, using the MCD from 
samples where the two active excitons are not properly 
resolved certainly leads to an underestimation of the cou- 
pling, since their contributions tend to cancel each other. 
In addition, we will see now that strong effects are ex- 
pected from electron-hole exchange within the exciton, 
when the giant Zeeman splitting is of the same order as 
the A-B splitting. 

B. s,p-d interaction in ZnO based DMS 
with electron-hole exchange 

We now include the effect of the electron-hole ex- 
change. We will see that the effective exchange inte- 
grals difference {No{a — (3))x, seen by the excitons can 
be much smaller than the difference of the exchange in- 
tegrals valid for free carriers, NQ{a — /3). We follow the 
description given by B. Gil et al. when analyzing the 
effect of the biaxial stress in epitaxial layers, and 
we add the giant Zeeman effect Ge and Gh- The posi- 
tion of the three excitons A, B, G, which can be optically 
active from symmetry considerations, is obtained by di- 
agonalizing a 3 x 3 matrix. This matrix mixes the two 
excitons active in polarization (one formed from the 
electron- hole states \s |) \p'^ 1), which gives rise to ex- 
citon A if A2 > 0, and \s t) i)), and the G exciton 
state \s i) t)- It reads (apart from an overall shift) s^S 

/-A2 7 \ 

7 A2 -V2A3 (12) 
\ ~V2A3 Al - 7 / 

where Ai is a trigonal crystal field taking into account 
the effect of the trigonal strain,— 7 is the electron-hole 
exchange energy. The giant Zeeman effect adds to that 

Gh-Ge \ 

Ge-Gh (13) 

Gh-Ge) 



The same matrix, with opposite giant Zeeman terms, 
applies in <7~ polarization. 

The most important effect of the electron-hole ex- 
change within the exciton, for the present study, is that 
it mixes the A and B states active in a given circular 
polarization. 

If spin-orbit coupling (A2) dominates over exchange 
(7), then one observes the A and B excitons formed 
on one of the electron-hole states, \s I) T)i and 
|s T) i) (respectively for A and B if A2 > 0). These 
are eigenstates of the c-axis spin projection operator, so 
that they are fully affected by the giant Zeeman effect : 
In that case, one would measure directly the spin-carrier 
couphng No{a — (3) on the plot of Fig. El In the opposite 
case, the two exciton states are formed on the symmetric 
and antisymmetric combinations of the previous states, 
resulting in a spin triplet (dark exciton) and a spin sin- 
glet (bright exciton), respectively. Both have vanishing 
spin components along the c-axis, so to the first order 
they do not exhibit the giant Zeeman effect in the Fara- 
day configuration. As visible in Fig. El this is not the 
case, but any intermediate configuration is possible. 

Fig. Iior a) displays the position of the A and B ex- 
citons, from Fig. Etb), as a function of the splitting of 
exciton A. It is clear that we observe an anticrossing be- 
tween the two excitons. It is reasonable to ascribe this 
anticrossing the the electron-hole interaction, and from 
the minimum distance, the value of 7 is a few meV. How- 
ever, the two excitons do not show the same curvature, 
which suggests that we cannot restrict ourself to the in- 
teraction between A and B only. 

Fig. llOf b) compares the experimental position (sym- 
bols) of excitons A and B, as a function of the applied 
field, and of exciton C, at zero field, to the calculated 
eigenvalues of the Hamiltonian [Eqs (|12|l and H13|l . solid 
lines]. The average spin projection of Co ions is calcu- 
lated using Eq A good fit is obtained for two sets of 
values of the parameters (which cannot be distinguished 
at the scale of the figure): 

Set 1: Al = 51 meV, A2 = 3.0 meV, A3 = 6.3 meV, 

7 = 3.4 meV, iVo(a - f3) = 0.8 eV, 

or 

Set 2: Al = 53 meV, A2 = -2.1 meV, A3 = 0.1 meV, 
7 = 3.1 meV, Noia ~ (3) = -0.8 eV. 
The value of Ai corresponds to a trigonal stress of 
13 kbar, according to Ref. f29l 

In both cases, we obtain a very good fit of the lowest 
exciton, and the fit for the second optically active exci- 
ton is reasonable, although not as good. The splittings of 
both excitons are almost opposite, as experimentally ob- 
served (note that if we take the parameters from Ref. l30l 
the splitting calculated for B almost vanishes). We ob- 
tain no information about the valence band ordering (a 
similar conclusion was claimed by Lambrecht et al^ in 
their analysis of the effect of strain) . This is due to the 
fact that we mainly observe an anticrossing between the 
two excitons A and B, which is a symmetrical process. 
Hence, our fit is mainly sensitive to the relative value of 
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FIG. 10: (color online) (a) Position of the A and B excitons, 
for the same sample as in Fig. I^j (with 0.4% Co), measured 
at different values of the applied field. The horizontal axis 
is the splitting measured on exciton A, plotted as positive if 
the (T^ line is at higher energy, (b) Position of the A and B 
excitons, as a function of the applied field, and of exciton C at 
zero field, for the same sample as in Fig. \^ (with 0.4% Co). 
Symbols are experimental data (blue for a~ , red for a^), lines 
are calculated (including for exciton A). 



{Ge — Gfi) and 7. We may simply notice that the hypoth- 
esis A2 > leads us to a much more isotropic spin-orbit 
coupling. 

The value obtained for the spin-carrier couphng, 
No\a — /9| ~ 0.8 eV, has to be considered as a first es- 
timate. A larger value of the electron-hole exchange 7 
reduces the sensitivity of the exciton to the spin split- 
ting, implying a larger spin-carrier coupling. It should 
be particularly interesting to measure the excitonic split- 
ting on samples with a higher Co content. 

This value is smaller than in many other DMSs. For ex- 
ample, in Cdi_a:Coa:Se, another wurtzite Co-based DMS, 
the exchange integral difference is 2.1 eVm^^ However, 
comparable or even smaller values of Ngla — P\ have been 



also reported in wurtzite DMSs, for example -0.26 eV in 
Cdi_2:Cra;S44 It is also in significant disagreement with 
theoretical predictions. Blinowski and co-workers^ have 
calculated Nq/3 = —3.2 eV in the case of Zni^^MnjjO. 
They have also pointed out that the absolute value of 
Nof3 is expected to be larger for Co than for Mn ions. 
Since Noa is usually about 0.25 eV and does not de- 
pend strongly on the host material,'*^ one could expect 
No{a — j3) > 3.4 eV in the case of Zni_3;Co2;0. A large 
exchange integral, No{a — (3) = 2.3 eV, has been pro- 
posed for Zni_a;Coa;0 nanocristallitesji^ but the details 
of the spectra supporting this conclusion were not given. 
If we assume a positive sign of NqP (which is not the usual 
case but is not excluded by theory, e.g., in Cdi-xCr^^S)^^ 
and Noa = 0.25 eV, then we obtain Nof] about 1 eV. 
Such a small value puts severe restrictions to the possible 
influence of carriers on ferromagnetism, at least through 
Zener-like mechanisms. 



C. Magnetism of Zni a;Co3:0 

A measure of the excitonic Zeeman splitting and of 
the MCD signal near the energy gap is known to repre- 
sent an efficient method to study the magnetic properties 
of DMSs. The dependence on the magnetization results 
from s,p-d interactions, which are strong with substitu- 
tional magnetic ions. These magneto-optical properties 
are less affected by the presence of secondary phases or 
inclusions, which fully contribute to the magnetization 
as measured by SQUID. Hence the observation of a clear 
effect of the magnetization on the excitonic or band-to- 
band magneto-optical properties of the semiconductor is 
a good hint, that the observed magnetization is due to 
the DMS. 

In the whole set of studied samples, both the Zeeman 
splitting and the near-bangap MCD signal agree with the 
idea of a paramagnetic system of isolated Co impurities 
(Maxwell-Boltzmann occupancy of the |±|) and |±i) 
spin states of Co^+ [see Eq. Q]. 

This model is based on the observation of character- 
istic d-d transitions of the 3d'^ electronic configuration 
of Co, which confirms the values of the parameters de- 
scribing the anisotropy and Zeeman effect in the ground 
state. The intensity of these lines increases linearly with 
the Co content, which suggests that all Co ions are in- 
corporated in substitution of Zn. A Maxwell-Boltzmann 
distribution over the spin components in the ground state 
well accounts for the intensity of the d-d transitions, and 
for the giant Zeeman effect of the optically active exci- 
tons. This demonstrates that Zni_j;Co2.0 is a DMS with 
a coupling between the bands of the semiconductor and 
the localized spins. 

It also stresses some differences with respect to the 
most usual case of Mn-based II- VI DMSs. The mag- 
netization of the single Co ion is not described by the 
isotropic Brillouin function used for cubic II- VI DMS 
such as Cdi_a;Mna;Te. The anisotropy causes a step in 
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the field dependence of the magnetization at low temper- 
ature, and a deviation from the Curie law in the thermal 
dependence of the low-field susceptibility. On samples 
with a few % Co, we did not need to introduce an effec- 
tive temperature T^f f — T + Tq in order to describe the 
influence of spin-spin interactions between neighboring 
spins (at least, beyond the nearest-neighbor spins). Any 
such phenomenological Curie- Weiss temperature Tq has 
to be smaller than 0.5 K for concentration up to 5%. 

On the other hand, the model assumes isolated Co ions, 
so that the eventual magnetic moments associated with 
Co pairs and complexes are significantly smaller than the 
moment associated with isolated Co ions. The negligi- 
ble contribution of anti-ferromagnetically coupled near- 
est neighbor pairs at low temperature can be explained 
by short range superexchange interaction integrals higher 
than a few tens of Kelvin. 

Contrary to the resonant peak close to the energy gap, 
the MCD signal related to internal Co transitions is not 
necessarily proportional to the Co magnetization. This 
signal is governed by the occupancy of the cobalt spin 
sublevels, and the Zeeman shifts and selection rules of the 
multiline absorption structure. Experimentally, its tem- 
perature dependence is also much weaker than the cor- 
responding one for the energy gap resonant MCD peak. 

In any case, we have not found any evidence of a fer- 
romagnetic behavior. Upon increasing the Co content, 
the excitonic Zeeman splitting or the MCD continue to 
follow the paramagnetic behavior expected for isolated 
Co impurities - with perhaps the exception of a small 
additional contribution below 5 T at our lowest temper- 
ature, which tends to increase with the Co content (com- 
pare Figs. [3 and EJ. No hysteretic behavior was detected. 
This contrasts with other observations on samples which 
have been grown by a different method, which exhibit 
ferromagnetic magnetization with an easy axis along the 
c-axis (an anisotropy which is opposite to the one we see 
for isolated, substitutional, Sd'* cobalt) ^ or with samples 
which have been annealedii^ However, it fully agrees with 
a study of the present (or similar) samples by magnetom- 
etry and electronic paramagnetic resonance.^® 

In a study of the magneto-optical properties of such 
samples, one should keep in mind the trend exemplified in 
Figs. [7| and (the presence of an additional contribution 
to the MCD signal below 5 T and 1.7 K), and also the 
strong effect of electron-hole exchange in the exciton. 



VI. CONCLUSIONS 

We observe absorption bands and lines due to Co^+ 
intra-ionic d-d transitions. The intensity of ^EE 
is proportional to the total Co concentration. We ob- 
serve reflectivity structures related to A, B and C exci- 
tons. Their energies increase with Co concentration. We 
measured the exchange Zeeman splittings of A and B 
excitons, which increase with Co concentration and de- 
crease with temperature. These exchange Zeeman split- 



tings and the intensity of energy gap resonant MCD sig- 
nal are described by using a model based on thermal 
occupation of |±|) and |±i) spin levels of purely para- 
magnetic Co^"'' ions. A first estimation of spin-exciton 
exchange integrals iVola — /?| was deduced, which results 
in a higher value of the spin-carrier integrals, of the or- 
der of 0.8 eV, when electron-hole exchange interaction is 
taken into account. 
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APPENDIX: REFLECTIVITY DUE TO 
EXCITONIC POLARITONS 

We present here the details of the polariton model used 
to fit reflectivity spectra. The model takes into account 
the strong coupling between excitons and photons in ZnO 
(the longitudinal transverse splitting, of several meV, 
is of the order of the excitonic linewidth). The dielec- 
tric function and reflectivity spectra near excitonic reso- 
nances are described in detail by Hopfield and Thomas^i 
and by Lagoisi^^i^ Only A and B excitons were taken 
into account. The dispersion of the polariton branches 
associated with a propagation of the light parallel to the 
c-axis was calculated numerically by solving of the fol- 
lowing equations (|A.1|I and \K.2\ : 



Lu\ — Lo'^ + [hk'^LjJA/m*) — ILOT A 

47raoBtJB 



(A.l) 



e(fc, Lo) — k c /lj 



(A.2) 



where eoo — 6.2 is the background dielectric constant of 
ZnOf^ aoA.B is the polarizability, F^ ^ are dumping pa- 
rameters, loa.b are the energies of excitonic transitions, 
and m* — O.STtoe is the effective mass of excitons."*^ 

We note ei, £2, £3 the three solutions of (|A.1|I and 
(|A.2|I . For each solution, we can determine the associated 
refractive index rii — y/ei, the wave vector ki = riiLojc 
and the polarization of A and B excitons: 



^A,B - + {fikfuJA,B/m*) - iujTA,B 



(A.3) 
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The reflection coefScient was derived from surface 
boundary conditions for electric field HA.4() . magnetic 
field I) A. 5(1 and the total polarization of A ((A.6p and B 
(|A.7|1 excitons: 

Ej+Er = E,+E2+E3, (A.4) 



Ei-Er^ niEi + n2E2 + ng-Eg, (A.5) 



With the derived expression for n\ we can calculate 
the reflection coefficient of the layer R = \ which 
takes into account the presence of the polariton branches. 
However, it is commonly assumed that the formation of 
excitons is not possible within a few nanometers close to 
the surface. A dead (exciton-free) layer with a thickness 
d (about two times the excitonic Bohr radius), and a 
background refractive index n = ^Jeq are usually added 
to the model. The reflectivity of the bilayer structure 
leads to an effective refractive index given by : 



'E2 



Eo 



4E3 - 0, 



'E3 = 0. 



(A.6) 



(A.7) 



{n^ + n)e^^'^ ~ n + v) 



(A.ll) 



Ej is the electric field amplitude of incoming light, Er is 
the corresponding one for reflected light, Ex, E2 and E^ 
are the electric field amplitudes of the three polaritons 
propagating inside the sample. From ((A.fi|l and (|A.7|I . 
we get: 



(A., 



and 



E2 


B A 
^1 ^3 


^1 £3 


El 


- A B 
^2 ^3 


^2 ^3 


E3 


B A 
^1 ^2 


^1 ^2 


El 


- A B 
^3 ^2 


^3 ^2 • 



Next, we obtain an effective refractive index 



Ei — Er 
Ei + Er 



ni 



"3t 



El 



E3 
El 



(A.9) 



(A.IO) 



Finally, the reflectivity spectra are given by : 



R = 



n* - 1 



71* + 1 



(A.12) 



In our samples, the reflectivity spectra are slightly 
modified by additional non-resonant absorption. We ob- 
tained a better fit of the experimental spectra by adding 
an imaginary part ie' to the background dielectric con- 
stant in Eq (|A.1|) . Examples of the fits of the re- 
fiectivity spectra of ZnO and Zni-ajCo^O epilayers are 
shown in Fig. |S1 
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